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To investigate the influence of orientational degrees of freedom onto the dynamics of 
molecular systems in its supercooled and glassy regime we have solved numerically the 
mode-coupling equations for hard ellipsoids of revolution. For a wide range of volume 
fractions (f) and aspect ratios Xq we find an orientational peak in the center of mass 
spectra x'oooili ^0 an d 0ooo(<?' Co O about one decade below a high frequency peak. This 
orientational peak is the counterpart of a peak appearing in the quadrupolar spectra 
Xwrnili^) an d 4>22m,i < l-i UJ )- The latter peak is almost insensitive on for xq close to one, 
i.e. for weak steric hindrance, and broadens strongly with increasing x . Deep in the glass 
we find an additional peak between the orientational and the high frequency peak. We 
have evidence that this intermediate peak is the result of a coupling between modes with 
/ = and 1 = 2, due to the nondiagonality of the static correlators. 

1. INTRODUCTION 

Although experimental data on plastic crystals [|I|-|3J have demonstrated that orienta- 
tional degrees of freedom by themselves show glasslike dynamics and even can undergo a 
glass transition their role for the dynamics and the glass formation of supercooled liquids 
has not been studied in great detail. With the extension of mode-coupling theory (MCT) 
from simple liquids to molecular systems P-II there is a microscopic theory which 



allows to investigate the role of the rotational degrees of freedom and their coupling to 
the translational ones. It is desirable to choose systems for which this coupling can be 
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varied systematically. From this point of view hard ellipsoids of revolution represent an 
interesting model system. Keeping their minor axis b fixed, a change of the aspect ratio 
Xq = a/b, where a is the major axis, will change the steric hindrance between the ellip- 
soids. This again will result in a change of the translational-rotational coupling. Applying 
MCT to hard ellipsoids it was shown that for xq about two or larger the ellipsoids freeze 
into a glass at a critical volume fraction <ft c {xo) due to the appearance of a medium range 



orientational order which originates from a precursor of a nematic order |12| . Recently we 
have shown that this precursor of nematic order also has an influence onto the dynamical 
features in the microscopic time and frequency regime [I3|. From both, the solution of the 
corresponding MCT-equations and from a MD-simulation a peak in the compressibility 
spectrum has been found at u op about one decade below a high-frequency peak (hf-peak). 
The additional peak clearly has an orientational origin since uj op scales with J -1 / 2 , / being 
the moment of inertia with respect to the minor axis. Therefore we called it an orienta- 
tional peak (op). These calculations were restricted to a single aspect ratio, x = 1.8, and 
one volume fraction, only. It is the motivation of the present contribution to study the 
microscopic dynamics and in particular the op as a function of xq and (ft — f &o%o where 
£0 is the number density. In this respect we would also like to mention that studying 
the athermal system of hard ellipsoids can also be instructive for real molecular liquids 
composed of linear molecules like e. g. diatomic molecules. Indeed, it has been demon- 
strated that the static correlators Si m j> m ' (q) and the non-ergodicity parameters jim.ym 1 (q) 
for diatomic molecules with Lennard- Jones interactions obtained from a MD-simulation 
n| can be well reproduced by a system of hard ellipsoids with xq = 1.5 [IS]. Note that 



this aspect ratio was not the result of fitting procedure, but could be determined from 
the Lennard- Jones parameters. 

Besides hard ellipsoids MCT has already been applied to study the dynamics of a 
dumbbell- molecule in an isotropic liquid of hard spheres |16[ , [[17]] and of a liquid of dipolar 
hard spheres [|TJ|]. Like for hard ellipsoids, the change of the elongation ( of the dumbbell 
allows to vary the steric hindrance with the hard sphere environment. This effect has been 
investigated quantitatively for the orientational correlators Ci(t) = (Pi(e(t) ■ e(0))) where 
Pi are the Legendre Polynomials and e(t) is the unit vector pointing along the symmetry 
axis of the dumbbell at time t. It has been found that weak and strong steric hindrance 
can lead to quite different relaxational behavior. For instance, in case of sufficiently small 
elongations the a-relaxation scaling law is strongly disturbed. The correlator with odd-/ 
and even-/ behave different. It is interesting that this finding agrees with that found from 
a MD-simulation of a diatomic Lennard- Jones liquid. JT9[ . 
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For dipolar hard spheres it is the dipolar interaction which is responsible for a translational- 
rotational-coupling. This type of coupling is different to that generated by steric hin- 
drance. Nevertheless, we have also found in the compressibility spectrum an op about 
a decade below a hf-peak |]IB| |. This may suggest that an op can be a rather general 



feature of a molecular liquid. This motivates to study for hard ellipsoids the microscopic 
dynamics as a function of the aspect ratio and volume fraction. 

The organization of our paper is as follows. In the next section we will introduce the 
basic quantities and their equations of motion. The results from a numerical solution of 
these equations will be presented and discussed in section III and in section IV we will 
summarize and will give some conclusions. 

2. Equations of motion: Hard ellipsoids 

We consider a system of N hard ellipsoids of revolution with mass M and moment of 
inertia / = M(l + Xq) /20 with respect to the minor axis, xq = a/b is the aspect ratio, i.e. 
the ratio of the major and minor axis. The length and time unit is chosen such that 6=1 
and M = 1. The microscopic configurations are given by the center of mass position x n 
and the set of Euler angles Q n = (<p n , $ n , Xn), n — 1, 2, . . . , N. Because of the rotational 
invariance around the symmetry axis of the ellipsoids of revolution the third Euler angle 
Xn is not involved in the steric hindrance. Consequently this degree of freedom will not 
be taken into account. Then, the basic quantities for the study of the molecular dynamics 
are the tensorial density correlators 

Sl m ,l'm'(q,t) = ■Jy(Pi m (q,t)pi'm'(q,0)) (1) 

where pi m (q, t) is the microscopic tensorial density determined by {x n (t)} and {{p n (t), $ n (t)}. 
I = 0, 1, 2, . . . and —l<m<l specify the rotational part of the density modes. For com- 
pleteness we mention that the dynamics of molecular liquids can also be described by a 
site-site formalism ||. However, hard ellipsoids are not made up for a finite number of 
sites. Hence one has to approximate a hard ellipsoid by, e. g. a two-site dumbbell where 
each site is the center of a sphere as it was used in refs. fT6fl , JE 



The correlation matrix S(q, t) = {Si m ji m '(q,t)) simplifies as follows. First of all we can 
choose the g-frame, i.e. it is q— (0,0, q) , q — \q\. This makes S(q, t) diagonal in m |]10||: 



Si m ,i'm'(q,t) — Su' m (q,t)5 mm i . (2) 
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In addition one can prove [JTOj that Su/ m (q,t) are real and that 



Su>- m (q, t) = S u > m (q, t) . (3) 

Second the head-tail symmetry of the ellipsoids implies that the collective correlators 
Sim,i'm'{(i,t) vanisch for I and/or I' odd. For its self part \ Vm , (q, t) , which will not be 
considered here, this is only true for I + V odd. 



How to derive equations of motion for S(q,t) has been described in great detail fLO 
Here we repeat just the result: 



S{q,t)+i £ cfS a (q,t)=0 (4a) 

a=T,R 



S a (q,t) +icfJ a ($S- 1 ((f)S(q,t) 
+3 a {q)jdt' m aa '(q,t-t')S a '(q,t') = (4b) 

a'=T,R 

with the initial conditions: 



S(£0) = S($), S a (g,0) = 0, 

S(g,0)=0, S«(?,0) = -zJ Q (g)q a (5) 

and 



(q a )zm,Z'm' — Sll'Smm' { r~T~ \ (^) 

y 1(1 + 1) a = R 



{J a {q))im,i'm' = {k B T/I a )8i V 5 mm > (7) 

where It — M and Ir = I. S(q) and J a (q) are the static density and current density 
correlators, respectively. The form of this set of equations is different but equivalent to 
that derived in ref . JTOj] . The main difference is the occurrence of the density - current 



density correlators S a (q, t) in Eqs. (f|). It has been demonstrated that the set of equations 
(£|) is appropriate for a numerical solution |18|]. Eqs. (£|) become closed by the MCT- 
approximation for the memory kernels "mfmi'm'iQi With the shorthand notation A = Im 
it is H 
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= ^EE (EK(?%i A "i^ 2 )c Al (?i) + (i - 2)])) . 

qiQ2 Aja; A" 
\ 2 \' 2 

•(^K'(gA'|g 1 A w ;g 2 A;)c v » Vi (g 1 ) + (1 <-> 2)]f S XlK {q x , t)S MX > 2 (q 2 , t) (8) 

A'" 

where u T (qX\qiXi, q 2 X 2 ) and u R (qX\q\Xi; q 2 X 2 ) are proportional to the Clebsh-Gordon 
coefficients C(lil 2 l; 000) and C(l\l 2 i, 101), respectively. They vanish for l + li + l 2 odd, and 
their explicit form is given in ||10|| . The important input into (§) is the direct correlation 
matrix c(q) = (c\x(q)) which is related to S(q) by the Ornstein-Zernike equation for 
linear molecules. 



c(g) = (47r/p )[l-S- 1 (g)] , (9) 
where po — N/V is the number density. 

Since we want to investigate the dependence of S(q, t) on the aspect ratio it is important 
to check whether the Eqs. (§) and Eq. (H) reduce for xq — 1 to the corresponding equations 
of motion for hard spheres. Due to the complete isotropy of hard spheres it is: 



Anc(q), (A,A') = (00,00) 
c\y(q) = S (10) 
^ =0, (A, A') ^ (00,00) 

where c(q) is the direct correlation function for hard spheres. Substituting Eq.(|H]) into 
Eq. (pi) strongly reduces the number of mode-coupling terms. First of all one obtains: 



m^,(q,t)=0, (a,a>)^(T,T) (11) 



Using u a (. . . ) from ref. 10 we get for the nonzero elements: 



/// 



rprp 
AA' 



/(&*) = Pl^^{(^-Qi c (Qi)) 2 S(qi,t)S X x'(q 2 ,t) + 



+ (-q- q 2 c(q 2 )) S X \>(qi, t)S(q 2 , t) + 



(12) 



+ \(q ■ Qi)(q- ft) c(qi)c(q 2 )[So 0i x'(qi,t)Sx,oo(Q2,t) + (!<-> 2)] J 
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where we used S(q,t) = Soo,oo(<f, t), the density correlator for hard spheres. Note that 
m^ T 00 (g, t) coincides with the memory kernel for simple liquids in general and for hard 
spheres in particular. But mjj,(q, t) seems to be nondiagonal in A. Such a nondiagonality 
would generate an unphysical influence (at least for hard spheres with smooth surfaces) 
of S\\(q, t) with A 7^ (0, 0) onto S(q, t). However, using the initial conditions Eq. @ and 



Sxx>(q) = , A^A' (13) 

which immediately follows from Eq. (P) and Eq. (|T0| ) one can prove that Eq. (|J) and 
Eq. (D yield: 



|^AA'(^)| t =o = , A^A' (14) 

for all v. Although Eqs. (^) and (§) are nonlinear we therefore expect that Eq. (|14|) 
implies 



S xy {q,t) = , A 7^ A' , (15) 

as for linear equations. In that case m^y (g, t) given by Eq. (|T2| ) becomes diagonal 
in A, i.e. the MCT-equations for the remaining correlators S\\(q,t), I > 0, decouple 
with respect to A. Therefore S(q,t) is not influenced by S\\(q,t) with A 7^ (0,0) and 
it obeys the MCT-equation for hard spheres. Because m AA (g, t) for A 7^ (0,0) contains 
bilinear terms S(qi,t)S\\(q 2 ,t) the hard sphere correlator S(q,t) has an influence on the 
orientational correlators S\\(q,t), which is unphysical, as well. The reason for this lies in 
our assumption that the density modes p\(q) are slow variables for all A, which only will 
be justified for x 7^ 1. Neglecting m(<f, t), Eqs. (|) reduce to a linear equation for the 
(normalized) correlators ®(q?t) = S~ 1/2 (g)S(g, t)S~ 1/2 (<f) @: 

*(q,t)+n 2 (q)*(q,t) = (16) 
with the hermitean frequency matrix squared: 



(ft (q))im, I'm' 



£( s ~ 1/2 (^) 



Im.V'm" 



l"m h 



k B T 



M 



q2 + MLi»(i» + r 



m" I'm 1 ■ 



(17) 
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These microscopic frequencies will play an essential role for our study of the microscopic 
dynamics. Since the static correlators Si mt i' m >(q) for hard ellipsoids are nondiagonal there 
will be a coupling between the various modes already on the linear level of Eq. flip]). This 
type of coupling has been absent for dipolar hard spheres due to the use of Wertheim's 
solution for c\y(q) for which the direct correlation functions are diagonal in I and /' |L8 



Therefore in contrast to our results for dipolar hard spheres we may expect additional 
features for hard ellipsiods. The microscopic frequencies will become " renormalized" close 



to the glass transition and particularly deep in the glass. In ref. [|T^] it has been shown 
that the microscopic frequencies in the glass follow from: 

with C%t% = linwm^'^t). 
3. Results 

To solve Eqs. (^) and Eqs. (|§p numerically one has to cut-off q and I. Since the numerical 
effort is considerable for hard ellipsoids we have chosen I < l co = 2 and 30 g-values which 
were distributed non-equidistantly between q min = 0.51 and q co = 40. 



More details about the numerical procedure can be found in ref. [18]. Furthermore we 
have restricted ourselves to xq > 1. Although no strict symmetry exists between xq > 1 
and Xq < 1 one may expect similar behavior for xq < 1. It would be interesting to check 
this point. The direct correlation functions needed as an input for the memory kernels 
were calculated from Percus-Yevick theory [^1]. The temperature enters into J a {q), only 



(cf. Eq. (0)). We choose ksT such that the thermal velocity (ksT / M) 1//2 equals one. 
This choice also sets the scale for the microscopic frequencies Qtm,Vm'<i) an d Qim,i>m>(q)- 
All results will be given for the g-frame. 

Before we come to the dynamical results we present in figure [I] the phase diagram. The 
liquid phase is below the solid line and bounded by the thick part of the dashed line. 
The remaining part corresponds to the glassy phase. This diagram contains two glass 
transition lines which cross at ((f), Xq) = (0.498,2.42). This crossing relates to a change 
of the mechanism for the glass transition. For 1 < xq < 2.42 it is the cage effect which 
mainly leads to a freezing into a glass whereas for xq > 2.42 a precursor of a nematic order 
becomes responsible for the glass transition []12"| |. Note, that this change at Xq = 2.42 does 



not happen suddenly but it is a crossover phenomenon. Therefore the numerical value 
x = 2.37 should not be over-interpreted. 



The correlators Sw m (q, t) for I < l co = 2 and q < q co = 40 were calculated from Eqs. (f|) 
and (§) for the points shown in figure |I[ From these results we have determined the 
correlation spectra (j} u , m (q,uj) and the susceptibility Xu'miQ^) — UJ( i>u l mi < li u )- Fig ure @ 
presents the compressibility spectrum for q = 4.2 and pairs of (<j), xq) on the liquid side, 
but close to the glass transition line (full symbols in figure [l]). Here we include xq = 1.8 
for which the static correlators were obtained from a MD-simulation We have chosen 
double-logarithmic scales on which the characteristic features become more prominent. 
For hard spheres, i.e. xq = 1, only one peak, which we called the high frequency peak 
(hf-peak), exists at Uhf for uo > 10 _1 . Now, making xq larger than one, we observe an 
additional peak at u op about one decade below the hf-peak. This peak is what we called 
the orientational peak (op) in the first section. It already appears for x = 1.1, the 
smallest value we have chosen besides Xq — 1. With increasing aspect ratio, up to our 
maximum value xq = 2, the position of the op almost does not change whereas Uhf first 
increases and then decreases. The most striking features are (i) the op becomes broader 
and (ii) the ratio of the height of the op and the hf-peak increases with increasing xq. 
This behavior can be understood as follows. Increasing xq enhances the steric hindrance 
which leads to a stronger translational-rotational-coupling. Therefore the orientational 
motion which is an almost free rotation with a single frequency for xq close to one couples 
stronger and stronger to the translational motion which enters into the compressibility. 
This leads to a more hindered rotational motion involving a broader band of rotational 
frequencies such that the op becomes broader and it leads to an enhancement of the 
orientational contribution at u op relative to the translational one at Uhf. Since cu^f does 
not change with / [13] it must be related to the translational motion. The fact that a 
distinct "fingerprint" of the orientational dynamics can be observed in the compressibility 
spectrum already for xq = 1.1 might be surprising. But figure || supplies an explanation. 
There we show x'u> m (4-2, oj) and <t>' u , m (4.2, u) on linear scales, for the volume fractions 
marked in figure [I] for x = 1.1. Because these quantities for (IV m) = (200), (221) and 
(222) strongly resemble that for (IV m) = (220) we only show data for (000) and (220). 
X220 (4.2, tu) and 2 2o (4.2, cu) exhibit one rather sharp peak at u> w 8. For Xq — 1 it 



follows from Eq. @ with Eq. (|T0|) that Si mj i> m r(q) = 8ui5 mm > for / > 0. Using this, q = 
and / = 1/10 (for x = 1) we get from Eq. flT7|) the frequency fi| ec = (10 ■ 2 ■ 3) 1 / 2 = 7.75 
for a free rotation of a sphere with I = V = 2. Therefore it is obvious that the peak in 
X220 an d 0220 & t ^'op ~ 8 comes from an almost free rotation. If the rotation is completely 
free then X22o( ( ?' u ) an d 022o(9> u ) are proportional to 5(u — il^ 00 ). Although the coupling 
of the rotational motion into Xooo ^ s ver y small for |1 — Xq\ <C 1 it is the almost 5-like 
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behavior of Xwoil^) an d 4>22oi ( li LO ) which can contribute significantly to Xoooili^) an d 
o ' oo (g, c<j). This contribution is more apparent for x'ooo than for O ' OO , as can be seen from 
figure [| Whereas the op at cu op pa fi 2 ree occurs in x'ooo f° r an studied volume fractions 
it only becomes present in O ' OO closer to the glass transition. The fact that uj op pa u' 
supports the rotational origin of the op. Two more features follow. First, the spectra for 
{IV m) = (000) depend sensitively on whereas those for {Win) = (220) almost do not 
change, and second uj op is practically constant in contrast to the position uihf of the hf-peak 
which moves towards higher frequencies with increasing 0. The latter observation implies 
that with increasing the glass becomes more stiff with respect to the translational than 
for the rotational dynamics, provided Xq is close to one. 

Now let us consider the results for the largest value for x we have studied. Figure £| 
presents the spectra (again on linear scales) for xo = 2 and the volume fractions indicated 
in figure [1]. The glass translation occurs at <f) c {xo = 2) = 0.545. Comparison of figure |] 
with figure |3| shows that all the spectra have become more structured and broader due to 
the increase of x . The rather narrow peak in x'220 (4.2, u;) and 22O (4.2, u) at to' op ~ 8 
for xo = 1.1 has almost disappeared. It occurs as a shoulder in x'220 an d as a peak in 22O 
deeper in the glass. Its shape and intensity depends on 0. A well pronounced peak also 
occurs in x'220 and 22O at u int between to pa 15 and uo pa 25, depending on 0. This peak 
is present in the liquid (0 < 0.54) and the glass (0 > 0.55). 

These two peaks have their counterparts in x'ooo an d 0ooo & t u op ~ an d uji nt ~ uj' int , 
although the peak in x'ooo at ^ ~ &int becomes more pronounced deeper in the glass, 
and has a shoulder-like behavior for OOO , only. Besides these two peaks a hf-peak exists 
above u ~ 25, for all volume fractions. We have checked the dependence of these peak 
positions on the moment of inertia I and have found that only uo op pa J op shifts with I. 
Therefore the peak at uo op pa uo' op) which is much less pronounced then for xq = 1.1 is the 
orientational peak. The origin of the intermediate peak will be elucidated below. 

It is also interesting to study the g-dependence of the spectra. This has been done for 
an intermediate aspect ratio xq — 1.5 and = 0.64, i.e. deep in the glass (cf. figure 0). 
The result is depicted in figure || for 0^ o (<7, w) and I = 0, 2. The g-dependence of Xiiq{o.i w ) 
is quite similar. Even more pronounced than for Xq = 2 we observe two peaks in 22O 
at uj' op pa 10 — 20 and at u) int ~ 35. Again, these two peaks have their counterpart in 
o ' oo (g,c<j) at Uop and uo int . In addition there is a hf-peak at oj^f- The dispersion of these 
peak positions depends on I. u op is practically g-independent whereas <J exhibits some 
g-dependence, and the position of the intermediate and the hf-peak shows dispersion for 
/ = and I = 2. 
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To understand the existence and the g-dependence of these peaks at least qualitatively 
we have calculated for xq = 1.5, <fi = 0.64 and I < l co = 2 the eigenfrequencies Q 3 (q) of 
the 2x2 matrix {£lio,i>o(q)), ^i(<?) = ^2i,2i(<?), ^2(9) = ^22,22(5) and the "renormalized" 
frequencies Q u (q), v = 1,2 and 3 which are the eigenfrequencies of the 3x3 matrix 
(fif "'o(<?)) (cf. pi). Since &? %(q) = for / and/or V equal to zero, (pffif vo (q)) is indeed 
a 3 x 3 matrix for I, /' < 2. The results for Q u (q) and £l u (q) are given in figure |]a and ^b, 
respectively. Figure |6|b also contains cu op (q), uj int (q) and u>hf(q) determined from 4>Q 00 (q, to) 
for xq = 1.5 and = 0.64, i.e. deep in the glass. Since Q\(q) cannot be distinguished 
from fl2{q) there are practically three branches in figure ^a. Because ^(q) — > for 
q — > it is the longitudinal acoustic branch, whereas the others are optic-like. These 
optic-like behavior also appears for the "renormalized" frequencies &i(q) and Cl 2 (q) which 
seem to reproduce at least qualitatively the results for oj op (q) and Ui nt (q). The same 
holds for Uhf(q) which is in phase with Cl^q) where &3(q) is the "renormalized" bare 
frequency Q^(q). Unfortunately, our results for Cl u (q) become less accurate for small q. 



For an expalantion see ref. [18]. Therefore we cannot definitely determine the dispersion 
for q — > 0. The intermediate peak only appears deep in the glass. In that case the 
nondiagonal static correlator S200 (?) becomes large. Therefore we believe that the origin 
of that peak comes from the additional coupling between the modes with / = and / = 2 
which already exists on the linear level (see discussion at the end of section 2). 

4. Summary and conclusions 

We have calculated the dynamics of a system of hard ellipsoids as a function of the 
volume fraction and the aspect ratio Xq. In the spectra for the center of mass dynamics 
we have found an orientational peak at u op , roughly one decade below a high-frequency 
peak. Due to its shift with the moment of inertia and its appearance in the spectra of 
the " quadrupolar" motions (quantities with 1 = 2) its origin is mainly of orientational 
nature. This peak becomes broader and its intensity relative to the intensity of the hf-peak 
increases with increasing xo, i.e. with the increase of the steric hindrance. Its position is 
practically q-independent and it scales with J -1 / 2 . This implies that it originates mainly 
from localized modes with q « 0, as already found in refs. and [|l8j. For the dipolar 



hard sphere system studied in ref. |18| it has also been speculated that this orientational 



peak could be one of several contributions to the boson peak. Since hard ellipsoids are an 
athermal system we can not draw such a conclusion here. However, the appearance of 
the op for all aspect ratios we have studied together with its existence for dipolar hard 
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spheres may give a hint that it may occur in many molecular liquids. That it is not an 
artefact of the approximate character of the used equations of motions has clearly been 



demonstrated by the MD-simulation [13 



It is also interesting for xq close to one that a variation of <fi almost does not affect the 
spectra for I — V — 2, in contrast to those for I — I' — 0. There is a significant influence 
of the rotational motion on x'ooo anc ^ ^ooo muc h less feedback from the translational 
motion to x'220 an d ^220 f° r wea k steric hindrance. 

Finally, deep in the glass we have found an intermediate peak between the op and 
hf-peak. This peak did not appear for dipolar hard spheres where the static correlators 
were taken diagonal in / and /'. Since this is not the case for the system of hard ellipsoids 
we believe that the intermediate peak turns up due to a coupling between modes with 
I = and I = 2 mediated by the nondiagonal elements of Sa> m (q). The q-dependence of 
the orientational, intermediate and the hf-peak in the compressibility spectrum can be 
described by the "renormalized" eigenfrequencies, at least on a qualitative level. 
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Figure 1. Glass transition phase diagram for hard ellipsoids: There are two glass transition 
lines (solid and dashed lines). See text for an explanation. The symbols show those 
points for which we have calculated the time-dependent correlators. To demonstrate the 
deviation from a hard sphere the cross section of an ellipsoid is shown for xq = 1.1 and 
x = 2.4. 
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Figure 2. Compressibility Xooo(4-2, for (0, xo) indicated by the full symbols in figure |l|. 
The result for hard spheres (xq = 1) with = 0.5184 is shown by the thick dashed line. 
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Figure 3. uo- dependence of xj' ffl (4.2,o;) (left) and </>/ ro (4.2, u) (right) for I = 0, 2, xq = 1.1 
and the 0-values marked in figure 1:0= 0.49 (solid) 0.50 (dotted) 0.51 (short dashed) 
0.52 (long dashed) 0.53 (chain) 0.54 (solid) 




Figure 4. cu-dependence of x^ (4.2,u;) (left) and <^ (4.2,u;) (right) for I = 0, 2 , xq = 2 
and the 0-values marked in figure ^ (ft = 0.53 (solid) 0.54 (dotted) 0.55 (short dashed) 
0.56 (long dashed) 0.57 (chain) 0.61 (solid) 
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Figure 5. w-dependence of (f)'{ l0 (q,Lj) for I = 0,2 ,x = 1.5,0 = 0.64 and q = 3.1 (solid 
line), 4.2 (dotted line), 4.7 (short dashed line), 6.5 (long dashed line) and 9.8 (chain line) 




Figure 6. g-dependence of (a) (solid line), fi^ (short dashed line), f2 x (long dashed 
line) and f2 2 (dotted line) (see text), (b) &i (dashed line), Cl 2 (dotted line) and f2 3 (solid 
line) (see text). The peak positions u op (diamonds), uj int (open circles) and u h f (full 
circles) were deduced from 0ooo('2S °°) f° r x o — 1-5 and = 0.64. 



